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ABSTRACT 


In the design of spacecraft for proper thermal balance, accurate information 
on the long-term optical behavior of the spacecraft outer skin materials is 
necessary. This paper presents a phenomenological model for such behavior, 
explains the underlying principles involved, and gives some examples of the 
model’s fit to actual measurements under simulated earth-orbit conditions. The 
paper concludes with some comments on applicability of the model to materials 
testing and thermal modelling. 


INTRODUCTION 


Since the early 1960 f s, a considerable body of data has been acquired on the 
performance of Alzak 1 anodized aluminum under simulated space radiation conditions. 
This material has generally poor space stability, but is unique in that it demon- 
strates a consistent, repeatable deterioration in optical properties when exposed 
to ultraviolet radiation in vacuum 2 . This characteristic makes it an ideal refer- 
ence material for validation of space simulation facility performance. 

In 1981, major modifications to the AFWAL/ML 3 SCEPTRE (Space Combined Effects 
Primary Test Research Equipment) began, with a goal of raising the then-inadequate 
(.2 UV sun) simulated solar intensity to a useful (>1.0 sun) level, and .1 mil 
Alzak was decided upon as a control material. By 1985, this effort had been con- 
cluded, which provided a unique opportunity to examine Alzak behavior under sim- 
ulated solar UV intensities in the range 0-1 sun. Examination of the historical 
record provided further information from an earlier work 4 which has extended this 
data set to 3 suns (SCEPTRE data at >1 sun intensity has been in a combined UV/ 
electron environment). Thus, data is available for an extremely wide intensity 
range under known conditions with a known material. 


^Research performed under Air Force Contract No. F33615-84-C-5130 , University of 
Dayton Research Institute 
*A product of Alcoa Corporation 

2 "Proceedings of the International Symposium on Solar Radiation Simulation" - 
IES/ASTM, January 1965, pps. 115-122. 

3 Air Force Wright Aeronautical Laboratories/Materials Laboratory 
4 "Effects of Vacuum-Ultraviolet Environment on Optical Properties of Bright 
Anodized aluminum Temperature Control Coatings," AFML-TR-67-421 , J. H. Weaver, 
AFML , May 1968. 
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The desire to understand the material behavior arises from the need for space 
simulation testing to be conducted under accelerated conditions, i.e. at inten- 
sities greater than those experienced in the actual flight environment. Testing 
is usually conducted at intensities for 2-10 suns, and then "scaled" by multiply- 
ing the test time against the intensity to yield what is presumed to be a useful 
dose measurement, generally called "Equivalent Sun-Hours" or ESH. The unstated 
assumption in such scaling is that of Reciprocity - that a reciprocal relation- 
ship exists between the intensity and the time constant for material damage, with 
no significant deviation in total damage to the materials 

The SCEPTRE tests of Alzak do not bear out this assumption. While data from 
individual test runs can be nicely fitted to a "french curve" asymptotic trace 
(Figure 1), when the data is expressed in terms of the total dose in ESH, any 
such relationship disappears (Figure 2) . This suggests that a more complicated 
relationship between time, intensity, and damage occurs. 

It is the intent of this paper to express, as completely as possible, that 
relationship. 


DERIVATION OF THE MATH MODEL 


Consider a monomolecular layer of Alzak containing some number N Q of possible 
sites for color center formation. Figure 3 illustrates the behavior of such a 
layer when exposed to ionizing radiation at an intensity of 3 ionizing photons/ 
unit area/unit time. One photon interacts with a site, forming a color center. 
Another photon interacts with an existing color center - and has no effect. The 
last photon passes through the layer without effect, presumably to create a color 
center deeper in the material. 


Our objective is to derive an equation which will define the concentration of 
color centers in this layer as a function of time. 

/N -N(t)' 




( 1 ) 


this equation is satisfied by N(t) defined as 

N(t) = N q (l-e _t/T ) ; t - ~ 

o 

which may be divided through by N to yield the concentration of color centers 
as a function of time 


( 2 ) 


N(t) _ . -t/x 

— - I-* 

o 


(3) 


The observable value related to the concentration of color centers is the solar 
absorptance (a g (t)), which may be found by 

a (t) = a (N -N(t))+a (N(t)) (4) 

s o o 1 


Subtraction of the initial value yields the change in solar absorptance as a 
function of time 
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(5) 


Aa (t) = a (t)-a 
s s o 


otj^d-e t//x ) 


Now this model is very interesting, despite its flaws, since this is the model 
implied by the assumption of linear reciprocity stated in the introduction. 

Indeed, a 100 hour exposure to radiation at 2 UV sun intensity should, according 
to this model, have exactly the same effect as 200 hours of exposure at 1 UV sun. 

Attempting to fit the observed data to this model yields Figure 4. Clearly 
the model fails - and the most obvious way in which it could fail is by assuming 
that all the damage occurs in a single monolayer of material. 

Therefore, let us consider a multi-layer model. 

Figure 5 illustrates such a model, which is essentially a stack of monolayers 
as defined in Figure 3, viewed from the side. Again, we see three incident pho- 
tons. One forms a color center in layer one, one impacts an existing center with- 
out effect, and the third passes through to create a new center in layer 2. The 
difference between this case and the monolayer case is now clear - at heightened 
intensities, damage will be done deeper in the material, so the total saturation 
value for color center concentration (and, by inference, solar absorptance) will 
be higher. 


The model is described by 




N -N(x,t) 


N 


I(x) ; I(x)=I Q e 


-ax 


which is satisfied by 


N(x,t) = N o (l-e' t/T(x) ) ; x(x) = 

Dividing through by yields the Color Center Concentration 

N(x,t) , -t/x(x) 

N 

o 


( 6 ) 


(7) 


( 8 ) 


This completes the derivation of the model, but does not yield an observable 
quantity. To get an observable (a g (t)) from this, we again apply Equation (4), 
noting this time that we must integrate over the thickness of the material, which 
leads to 


f 


a (t) = a I [N -N(x,t) ]dx+a 


r D 

1 N(x, 

1 Jo 


t)dx 


f 


(9) 


N dx 
o 

rearranging terms, and integrating the simpler ones yields 

* D 


a (t) = a N D+(a 
s v o o 


rV f\ 


. -atl e ax )dx 
(1-e o 


N D 
o 


( 10 ) 
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This cannot be analytically solved for all time, but it can be solved for the 
trivial case of t=0 (yielding a , as expected), and for the limit as time goes 
to infinity. This last yields 


lim 

t~*oo 


[a s (t>] 


a N D+(a -a n )N d 
o o 1 0 o 1 

N D 
o 


(ID 


Where we note that at depths greater than d^, no color centers are formed, so the 
integral range for the N(x,t) term is thus restricted. Cancelling the N ! s 

« 1 j o 

yields 


lim 

t-x» 


[« s (t)] 


“o +( °r“o ) t 


( 12 ) 


And subtraction of the yields 


Aa s = a s' a 0 = (a r a 0 ) ~D 


(13) 


Which expresses the total change in absorptance over infinite time. It merely 
remains to find the maximum depth of damage d . We do this by applying Beer’s Law 
to compute the intensity at the maximum depth. 


Taking the log of both 
maximum damage must be 
ically improbable that 


I(x) 


x=d 


1 


I(cL) = I e ad l 
1 o 


(14) 


sides, noting that the threshold intensity at the depth of 
a constant (since at greater depths it is quant um-mechan- 
there will be any effect) and solving for d^ yields 


d = A In (I ) + B 
1 o 


(15) 


where A and B are arbitrary constants. Substituting this back into Equation (12) 
and combining constants yields 


^ Act s (t) = (a 1 -a Q )(A ln(I Q )+B)/D (16) 

= A In (I )+B 
o 

We now make the assumption that the complete, time-dependent equation can be 
separated into static and dynamic components, of which Equation (16) is the 
static component 


Aa g (t) « [Aa(t)]f(t) (17) 

t=°° 

= (A ln(I o )+B)f(t) 

Inspection of Equation (10) yields the dynamic (time-dependent) component 

f(t) = K f X [l-e -t/T(x) ]dx; t ( x) = at o e~ ax (18) 

^ o 
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Which has no analytical solution. We approximate this function by converting it 
to series form 


n 


f(t) = K £ (1-e 

i=l 


-t/T (x) 


)Ax ; x=Ax(i-l); 

/ n T -aAx(i-l) 
t(x)=gI e v 
o 


(19) 


And then note that the most significant term is the first, due to the nature of 
the exponential component. Truncation of all terms for i>l yields 


f(t) = K(l-e t/x ) ; t - (20) 

o 

where K is an arbitrary constant. Substitution back into Equation (17), and 
combination of the constants finally yields 

Aa s (t) = [Aa s (t)] oo f(t) = (I q ) (l-e" t/T (I o ) ) 


T(I ) = 
o 


1 

al 


_k_ 

I 

o 


a- (I ) - A In (I )+B 
1 o o 

where a, A,B = constant (21) 

Which represents an approximate expression for the observable variable - the 
change in solar absorptance. 


TESTING THE MODEL AGAINST REALITY 


Figures 6 and 7 illustrate that, indeed, t is inversely proportional to 
intensity, and a^ is linear with the log of the intensity. Simple least-squares 
fitting produces the constants in question, thus 

T = 95.0/1 

a = .008 log(I) 4- .042 

(or k=95, A=.008, B=.042 in the nomenclature of Eq. (21) 

Substituting these into Equation (21) and plotting for the actual test intensities 

of .2, .5, 1.0, and 3.0 suns yields Figure 8 - which shows that the equation does 

indeed fit the data to within the illustrated error limit of ±.005 in a . 

s 

This data fit is far too close to be fortuitous - if this model is not 
accurate, it is at least much closer than the naive model illustrated in 
Figure 4. 
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CONCLUSION - IMPLICATIONS ARISING FROM THE MODEL 


Given that Equation (21) provides a working approximation to material 
behavior in this case, we are left with several interesting questions: 

1 - Does the model apply to other materials? 

We do not know. Existing SCEPTRE data sets for materials other than Alzak 
encompass at best two exposure intensities, which is insufficient to validate 
linearity of a^ and x. However, if an iterative solution for a^ and x is toler- 
ated, every data set we have , for organic and inorganic paints , OSR* s , etc . can 
be fitted to Equation (8) . This procedure yields a resultant function which is 
valid only at the intensity for which data is available, but it is nonetheless 
encouraging. 

2 - Does the model apply to other radiation environments? 

Again, insufficient data prevents us from making a definitive answer, but 
the iterative procedure mentioned above has so far been successful at producing 
a fit to within experimental error limits for all SCEPTRE UV/vacuum and UV/vacuum/ 
electron tests. 

3 - Assuming that the model is (at least approximately) valid, what does it 

say about present-day accelerated testing of spacecraft materials for 
space stability? 

The model has several implications for accelerated testing, well expressed 
by comparing Figures 4 and 8. The assumption that test results can be corrected 
for varying intensity by scaling the time axis is clearly false . It does appear, 
however, that testing at accelerated intensities represents a "safe worst case”, 
since the saturation value is always higher for higher intensity. 

Perhaps the most important implications are for intercomparison of ground- 
based simulation data to flight test results. When those of us in the simulation 
community speak of a "UV-sun" we mean an equivalent UV intensity to that from the 
sun , assuming a view factor of 1 . 0 with the specimen normal to the incident 
light ray . In the real environment, for most any condition, view factor consider- 
ations will lead to a time-averaged intensity of less than .5 UV-suns for a real 
spacecraft on-orbit. If the simulation testing has been conducted at 2 suns, then 
the acceleration factor will be 4, not 2 as is often claimed. Given the intensity 
dependent nature of Equation (21), this suggests that such comparisons need to be 
examined very carefully - and leads one to question exactly what is meant if ground 
test results at high acceleration factors do indeed match flight test data. Given 
the prediction that tests at accelerated intensity will always yield more damage 
than a lower intensity, one is left to wonder which non-modeled component of the 
flight environment is responsible for the implied additional damage, and which 
materials are susceptible to its effects. 
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If Equation (21) does prove applicable to other materials and other environ- 
ments, then it may become possible to answer some of these questions, and to put 
modern ground-based space environmental effects testing on a more sound theoretical 
footing. In this regard, it should be noted that the models defined by Equations 
(1) - (21) are in fact special cases of a model defined by Equation (22) , 


3 N(x,t) 

at n 

o 


a(x,t) 


N -N(x,t) 

K*.t) 

O 


( 22 ) 


In which both the incident radiation intensity and the cross-section are defined 
as functions of both depth and time. If solutions - whether exact or approximate - 
can be found for this model, then it may well prove possible to generalize the work 
described herein to more complex systems and conditions. 
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SYMBOLS 


t 

N(t ) 

N 

o 

N(t) 

total 

oc 

S 

oc 

0 

cc 

1 

a 

k, A,B,a^ 

T 

I 

0 

l. 

1 

n 


i 

Ad 



a 


D 


time (hrs) 

number of cells in a given monolayer changed up to time=t 

number of cells per monolayer 

total number of cells changed at time=t 

solar absorptance (absorptance units) 
initial solar absorptance (absorptance units) 
final solar absorptance (absorptance units) 
radiation effectiveness cross-section (fraction) 
arbitrary constants 
time constant (hrs) 
surface radiation intensity 
radiation intensity at layer i 

number of layers to which radiation penetrates 
index of a given layer 

thickness of a single layer of material 
maximum depth to which radiation penetrates 
optical density of material at radiation wavelength 
total thickness of material 
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Figure 3. Degradation phenomena in one monolayer 
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Figure 5. Color center generation in depth 


205 



206 


Figure 6. Inverse linearity of t with intensity 
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Figure 8. Observed data fit to math model 



